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Abstract 



We propose an improved iterative scheme for calculating higher genus contributions to 
the multi-loop (or multi-point) correlators and the partition function of the hermitian 
one matrix model. We present explicit results up to genus two. We develop a version 
which gives directly the result in the double scaling limit and present explicit results 
up to genus four. Using the latter version we prove that the hermitian and the complex 
matrix model are equivalent in the double scaling limit and that in this limit they are 
both equivalent to the Kontsevich model. We discuss how our results away from the 
double scaling limit are related to the structure of moduli space. 
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1 Introduction 



The hermitian one matrix model is interesting for many reasons. The matrices involved 
being N x N, the model can be solved in the limit when — > oo, as was first demon- 
strated in Ref. Furthermore the model admits a 1/A^^-expansion. This expansion is 
the simplest example of the 't Hooft topological (or genus) expansion — in the given 
case of a D = dimensional quantum field theory. If the leading order (planar or genus 
zero approximation) result for some quantity, say the vacuum energy, is normalized to 
be of order 1, the genus g contribution is of order 0{N~'^^). Even though the hermitian 
matrix model is often considered a toy model, its perturbative expansion reveals the 
combinatorial problems one will encounter in the multi-dimensional case: The sum of 
graphs contributing to a given order of the genus expansion is convergent, while the 
attempt of summing over all genera leads to a divergence. This was one of the original 
motivations for studying the genus expansion of the hermitian matrix model. 

One more application of the hermitian matrix model is due to the fact that it 
describes discretized 2-dimensional random surfaces and hence 2D quantum gravity. 
In this language a diagram of genus g is equivalent to a (piecewise linear) surface of 
genus g. The hermitian one matrix model was extensively studied from this point of 
view after the observation that it might not only describe pure 2D gravity but also 
some minimal conformal field theories coupled to 2D quantum gravityi. The relation 
given in |^ for genus zero was extended to any order of the genus expansion in the so- 
called double scaling limit, where the coupling constants approach some critical values 
in a iV-dependent way [|]. 

In the double scaling limit many interesting approaches have been developed in 
order to deal with higher genus contributions. The technique of exactly integrable 
systems (KdV and KP hierarchies) was applied to the model in Ref. . It provides an 
algorithm for genus by genus calculations in the double scaling limit. Another approach 
proposed by Witten |0] is based on the interpretation of the double scaling limit of the 
matrix model as a topological theory. The problem is then reduced to calculating 
intersection indices on moduli space. In this approach it was possible to obtain explicit 
results in the double scaling limit up to genus four |Q. A further development of 
this approach was the representation of the generating function for the intersection 
indices as a hermitian one matrix model in an external field (the Kontsevich model). 

For the hermitian matrix model away from the double scaling limit, the original 
genus zero solution was obtained by solving the large-N integral saddle point equa- 
tion for the spectral density. To extend the results to higher genera, the orthogonal 

^Although the relation conjectured in [Q turned out to be wrong, the paper was nevertheless 
instrumental for the later developments. 
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polynomial technique was applied in Ref. |TD[. While this method works well for the 
partition function in the case of simple potentials and lower genera, it is difficult to 
implement if one wants an explicit calculation of multi-point correlators and if the 
potential is "generic". 

Another method of solving the hermitian matrix model is based on the loop equa- 
tions [0 which were first written down for the hermitian matrix model in Ref. [Q. In 
this approach one solves the genus zero equation and uses the solution in an iterative 
genus by genus procedure. In the scheme originally proposed in Ref. and elaborated 
in Ref. |jl3|, one determines the higher genera contributions by solving algebraically the 
higher loop equations and imposing on all correlators the same analyticity structure 
as that of their genus zero versions. Following this line of action genus one correlators 
for (j)'^ and (p^ potentials were calculated in ref. [14|. However, this method involves the 
entire chain of loop equations and is in practise only applicable when one only wants 
to calculate a few terms in the genus expansion for a potential with a small number of 
terms. 

An alternative iterative procedure for calculating higher genus contributions to the 
correlators of the the hermitian matrix model has been proposed recently fl^. It is 
based only on the first in the chain of loop equations, is applicable for an arbitrary 
potential and can be pursued order by order of the genus expansion without any prob- 
lems. A key point in this scheme is a change of variables to the so called moments 
of the potential. The scheme is close in spirit to the iterative solution of the unitary 
and hermitian (with cubic potential) one matrix models in an external field |T^i. The 
genus one partition function and correlators were calculated for an even potential in 
Ref. |jl8| and for an arbitrary potential in Ref. [^, where it was also proven that genus 
g contribution to the n-loop correlator depends on at most 2 x [3g — 2 -|- s) lower 
moments. An analogous algorithm for the complex matrix model was developed in 
Ref. |jl9| and the genus 2 and 3 correlators were explicitly calculated for an arbitrary 
potential. 

In the present paper we propose an improved algorithm for an iterative calculation 
of higher genus contributions to the multi-loop correlators and to the partition function 
of the hermitian one matrix model. We perform calculations up to and including genus 
two. Our algorithm differs from that of Ref. [|15| by a redefinition of the moments, 
which has the advantage that it allows us to develop a version which gives directly 
the result in the double scaling limit where we perform explicit calculations up to and 
including genus four. We compare the analysis to that of the complex matrix model 
and prove in particular that the hermitian and complex matrix models are equivalent in 



^For the Kontscvich model similar moments has been used in Ref. fVi 
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the double scaling limit. We formulate a limiting procedure which allows us to obtain 
the Kontsevich from the hermitian matrix model and use it to prove in a very direct 
manner the equivalence between the double scaling limit of the hermitian matrix model 
and the Kontsevich model. 

We discuss a possible interpretation of our results for the hermitian matrix model as 
representing a particular discretization of moduli space which preserves basic geometri- 
cal properties such as intersection indices, and which allows for a detailed study of the 



boundary of moduli space ||20[. We conjecture that the moment-technique developed 
in the present paper might be an efficient tool in the study of this "fine structure" of 
moduli space. 

The paper is organized as follows: In Section 2 we introduce the basic ingredients 
of the iterative procedure, the loop insertion operator, the loop equation and the new 
variables — the moments. In Section 3 we describe the iterative procedure and present 
the results for the average of loop operators and the partition function for the hermi- 
tian matrix model up to genus two. In Section 4 we demonstrate in detail how our 
iterative procedure works in the double scaling limit and calculate the correlators and 
the partition function up to genus four. In Section 5 we show that the Kontsevich 
model can be obtained as a certain limit of the hermitian matrix model, prove the 
equivalence of the double scaling limit of the hermitian matrix model and the Kontse- 
vich model and discuss the possible connection between our results for the hermitian 
matrix model away from the double scaling limit and the structure of the discretized 
moduli space. Finally Section 6 contains a short discussion of future perspectives and 
in the Appendix we develop a double scaling version of the iterative procedure for the 



complex matrix model proposed in ||l9[, and prove that the double scaling limits of the 



complex and hermitian models coincide. 



2 Main definitions 

The iterative procedure we are going to describe is based on three main ingredients, 
the loop insertion operator, the first in the chain of loop equations of our model and 
a suitable change of variables. It allows us to calculate the genus g contribution to 
the n-loop correlator for (in principle) any g and any n and (also in practise) for any 
potential. The possibility of going to an arbitrarily high genus is provided by the loop 
equation while the possibility of considering an arbitrarily high number of loops is 
provided by the loop insertion operator. The change of variables allows us to consider 
an arbitrary potential. 
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2.1 Loop insertion operator 

The hermitian one matrix model is defined by the partition function 

Z[g,N]=e^'^ = f d(j)exp{-NTTV{(j))) (2.1) 

where the integration is over hermitian N x N matrices and 

°° o- 

V{<P) = Y.-<t^ (2-2) 

j=i J 

Expectation values (or averages) are defined in the usual way as 

{Qm = \j rf0exp(-iVTr 1^(0)) Q(0) . (2.3) 
We introduce the generating functional (the 1-loop average) 

1 oo 
fc=0 

and the s-loop correlator (s > 2) 

oo 

iy(pi,...,p,) = iV'^-2 ^ (Tr/^..Tr0'=^)„/p^+^..p^+i (2.5) 

k\,...,ks=l 

where conn, refers to the connected part. One can rewrite the last two equations as 
follows 

W{p,, ...,Ps) = N'-' ( Tr . . . Tr . (2.6) 

As is explained in Ref. [|], these quantities are associated with the (Laplace images of 
the) sum over discretized open surfaces with s boundaries. 

The correlators can be obtained from the free energy, F, by application of the loop 
insertion operator • 



dV{p) ■ 



Ti^/ N d d dF 

where 

dV{p) p^+i dg, • ^ • ^ 

Notice that Eq. ( |2.7| ) can be rewritten as 
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which shows that if the 1-loop correlator is known for an arbitrary potential, all multi- 
loop correlators can be calculated. This is why it is named above as the generating 
functional. 

With the normalization chosen above, the genus expansion of the correlators reads 

OO 

Wip„ ...,p^) = Y^ W,{pu ...,ps) is>l). (2.10) 

9=0 

Similarly we have 

OO -| 

3=0 

for the genus expansion of the free energy. 
2.2 The loop equation 

The first in the chain of loop equations for the hermitian matrix model can conveniently 
be written as IPII 



r ^ YMwiu) = {W{p)r + -^,W{p,p) (2.12) 

where V{uj) = J^jOj^'^/j and C is a curve which encloses the singularities of W{u) 
and not the point uj = p. This contour integration acts as a projector picking up the 
coefficients of p^^. Due to Eq. (|2.9|) the second term on the right hand side of the 
loop equation ( |2.12|) is expressed via W{p), so that ( |2.12| ) is a closed equation which 
determines this quantity. 

To leading order in one can ignore the last term in Eq. ( |2.12| ) and assuming 

that the singularities of W{uj) consist of only one cut on the real axis and that W{p) 
behaves as 1/p as p oo one finds |jlT| 

2 Jc zm p — UJ [[UJ — x)[uj — y) j 

where x and y are determined by the matrix model potential in the following way 

do; yV^^O^ fdu^ -Vjuj) 

c27rt ^ (^^ _ x){uj - y) Jc 2tii 

This genus zero solution will be used below in the iterative procedure of solving the 
loop equation. 

Inserting the genus expansion (p.lO|) in Eq. ( p.l2| ) it appears that Wgij)) , g > 1, 
obeys the following equation 

{k - 2Woip)} W,ip) = E WAp) W.^Ap) + ^^^^'-i(P) (2.15) 
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where i^' is a linear operator, namely 

Knp)^i^rMf(.). (2.16) 

Jc 2m p — uj 

In Eq. ( p.l5| ) Wg{p) is expressed entirely in terms of Wg.{p), Qi < 9- This is what 
makes it possible to develop the iterative procedure mentioned in the Introduction. 

2.3 The new variables 

To characterize the matrix model potential we introduce instead of the couplings gj 
the moments and Jk defined by 

= f^l ^ W5 

These moments depend on the coupling constants (^^'s both explicitly and via x and y 
which are determined by Eq. ( |2.14| ): 

Mk = gk+ {^x + {k + l)y)gk+i + ... , (2.19) 
Jk = gk+ {{k + i)x + ly) gk+i + ... . (2.20) 

Notice that Mk and Jk depend exphcitly only on gj with j > k. 

There are several motivations for introducing these new variables. First, as we 
shall see below, for each term in the genus expansion of the partition function and the 
correlators, the dependence on the infinite set of coupling constants will arrange into 
a simple function of a finte number of the moments. Moreover, these new variables 
reflect more directly than the coupling constants the possible critical behaviour of the 
matrix model. Let us briefly describe how this comes about. 

Performing the contour integral in (|2.13| ) by taking residuals at uj = p and uj = oo 
one finds 

Woip) = \ [V'{p) - M{p)^{p-x){p-y)^ (2.21) 

where M{p) is a polynomial in p of degree two less than V{p). As already discussed in 
the Introduction, Wq{p) can also be determined by an analysis of the matrix model in 
the eigenvalue picture |l[]. Requiring that Wq{p) is analytic in the complex plane except 
for a branch cut [y^x] on the real axis and behaves as 1/p for p ^ oo corresponds to 
requiring that the eigenvalue density, p(A), has support only on the interval [y, x] and 
is normalized to 1. The eigenvalue density is in this situation given by 

p{X) = ^M{X)^{\-y){x-\) y<X<x. (2.22) 

71 * 
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This function vanishes under normal circumstances as a square root at both ends 
of its support. Critical behaviour arises when some of the roots of M(A) approach x 
or y. For a non-symmetric potential the so called m*^ multi-critical point is reached 
when (m — 1) extra zeros accumulate at either x or y. Comparing ( |2.13D and ( [2.21| ) it 
appears that 



fc-i 



Jk oc 



k-l 



(2.23) 



\=y 



so the condition for being at the m multi-critical point is simply Mi = M2 = . . . = 
Mk-i = and Mk 7^ if the extra zeros accumulate at x and Ji = J2 = ■ ■ ■ = Jk-i = 0, 
Jfc 7^ if the extra zeros accumulate at y. For the symmetric potential {V{(f)) = V{—(f))) 
the eigenvalues will be distributed symmetrically around zero and hence confined to 
an interval of the type [— ^^e one arc case. In a situation like this the m^^ 
multi-critical point is characterized by the eigenvalue density having (m — 1) extra 
zeros at both —^/z and y/z. Reexamining ( p.l7|) and ( p.l8| ) it appears that for the 
symmetric potential 



(2.24) 



Hence in this case we have only one set of moments and the condition for being at the 



m 



th 



multi-critical point is the vanishing of the first (m — 1) of these. 



This formalism obviously allows for a treatment of a more general situation where 
m extra zeros accumulate at x and n extra zeros accumulate at y. Such multi-critical 



models have been studied in Ref. |22|. We will however restrict ourselves to the tradi- 
tional models. 



The superiority of the moments defined in ( p.l7| ) and (|2.18|) as compared to the 
coupling constants will become even more clear when we consider the double scaling 
limit in Section |4.1| . By then it will also become evident why these new moments are 
more convenient than those originally introduced in Ref. [|15| . 



3 The iterative procedure 

Our iterative solution of the loop equation results in a certain representation of the free 
energy and the correlators in terms of the moments. In this section we describe the 
structure of Wg{p) and Fg and show that the iterative procedure can be conveniently 
formulated by referring to the eigenvectors of a linear operator, K. We prove that 
the genus g contribution to the s-loop correlator depends at most on 2 x (3g — 2 + s) 
lower moments (2 x {3g — 2) for the partition function) for g > 0. We perform explicit 
calculations up to genus two. 
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3.1 The structure of Fg and Wg{p) 

Our iterative procedure of solving the loop equation results in the following represen- 
tation of the genus g contribution to the free energy 

= E («i /3i ■■ • Ala, 7). ^"^ • -^J!";/;; • • • 9>i, (3.1) 



/3i>l 



where d = x — y is the distance between the endpoints of the support of the eigenvalue 
distribution, the brackets denote rational numbers and a, j3 and 7 are non- negative 
integers. The indices cti, . . . , a^, . . . , take values in the interval [2, — 2] and the 
summation is over sets of indices. In particular Fg depends on at most 2 x ["ig — 2) 
moments. Furthermore, since nothing allows us to distinguish between x and Fg 
must be invariant under the interchange of the two. Hence one gets 



F : {ai . . . Us] (3i . . . 7) = (-l)'^(A • • • A; «i • • • «, 7) • 



(3.2) 



There are certain restrictions on the integers which enter Eq. (3J.). Let us denote 
by Nm and Nj the total powers of M's and J's, respectively, i. e. 



Then it holds that 



and 



Nm = s-a, Nj = 1- p. 



Nm < 0, Nj<0 



(3.3) 
(3.4) 
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NM + Nj = 2-2g, (3.5) 
Fg : E-=i(«. - 1) + - 1) + 7 = 4(7 - 4 . (3.6) 

The relation (|3.5| ) follows from the fact that the partition function Z = e^o " 



IS in- 



variant under simultaneous rescalings of and the eigenvalue density, p(A); N ^ k ■ N, 
p(A) —>■ ■^p(A). The relation (|3.6|) follows from the invariance of Z under rescalings of 
the type N k'^ ■ N, gj k^~'^gj. Finally the following inequality must be fulfilled: 



s 
i=l 



l) + E(/5.-l)<3^7-3 



(3.7) 



This requirement becomes more transparent when we consider the double scaling limit 
in Section [4.1| . In combination with Eq. ( |3.6| ) it gives 



g-l<-f<4g-4. 



{3.t 
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To explain the structure of Wg{p), let us introduce the basis vectors x^"'\p) 



X^-Hp) = -^\^^:'\p)-j:X^'\p)Mn-k+i\ n>l, (3.9) 



k=l 

n-1 



¥^\p) = 1 |$W(p)-^^M/W(p)J„_,+i| n>l (3.10) 



where 



$W(p) = [p - {{p - x){p - y)}-^^^ n>0, (3.11) 

^^;\p) = ip~yr' {iP~x)ip-y)r'^' n>o. (3.12) 

It is easy to show for the operator K defined by Eq. (|2.16|) that 

{k-2Wo{p)}x^''Hp) = ^>1' (3-13) 

{K-2Woip)}¥-\p) = ^— i-^ n>l (3.14) 

and that the kernel of [k - 2Woip)} is spanned by $i°)(p) = 

Since Wg{p) can be obtained from Fg according to Eq. (|2.7]), the representation 
( pAl) implies§ 

Wg{p) = E + D^;^¥^\p)} g>l. (3.15) 

?i=i 

We do not add any multiple of $^°^(p) or ^^\p). Doing so would contradict the bound- 
ary condition W{p) 1/p for p ^ oo since this behaviour was already obtained for 
genus zero. We note that this structure of Wg{p) is in agreement with the assump- 
tion that Wg{p) is analytic in the complex plane except for a branch cut [y,x] on 
the real axis. 

The coefficients A^J^^ are of the same structure as Fg and the relation ( p.4| ) still 
holds. However in this case the indices ai, . . . ,as, Pi, . . . , Pi take values in the interval 
[2, 3g—n]. Hence Wg{p) depends on at most 2x [3g — l) moments. The invariance of the 
partition function under the rescalings described above has the following implications 
for the structure of A^J^^: 

4"^: NM + Nj = 2-2g , (3.16) 

: E-=i(«. - 1) + E;=i(/5, - 1) + 7 = % - 2 - n . (3.17) 



similar formula is proven in Ref. for a different definition of the moments and of the basis 
vectors. 
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We also have an analogue of (|3.7[) for v4^"\ It reads 



k I 

4") : ^(«, - 1) + 5] - 1) < 3(7 - n - 1 . 

i=i j=i 



(3.18) 



Again we refer to Section |4.1| for further explanations. However, we note that combin- 
ing (|3.17|) and ( |3.18|) one gets again the bound ( p.8|) on 7. 

As was the case for Fg, Wg{p) must be invariant under the interchange of x and y. 
This means that Df^ must appear from A^^^ by the replacements d — > — (i, J ^ M. 
(We note that we do not have a relation like ( ^.2| ) for the y4^"'''s.) Furthermore, Wg{p) 
must be an odd function of p for a symmetric potential. This implies that A^") = 

That the structure of Fg and ^^(p) actually is as described in this section can 
be proven by induction. We will not go through the proof here. Instead we refer to 
Refs. fl^, |19| . In the latter reference a formula somewhat similar to (|3.15|) was proven 
for the complex matrix model. However the strategy of the proof will be evident from 



Sections |3]^ and ^]3| where we describe the iterative procedures which allow us to 
calculate Wg{j)) and Fg for any g starting from Wq{j)). 

3.2 The iterative procedure for determining Wg{p) 

According to Eq. ( p.l2|) , we need to calculate Wq{p^p) in order to start the iterative 
procedure. To do this we write the loop insertion operator as 



d 



d 



dx d 



dy d 



dV{p) dV{p) dV{p)dx dV{p)dy 



where 



d 



dV{p) 



E 



2 ^ 

1 p'+^ % 



The derivatives dx/dV{p) and dy/dV{p) can be obtained from ( |2.14] ) and read 



dx 



Using the relation 



dV{p) Ml 
d 



dy 



one finds 23 



dV{p) 

Wq{p,p) 



dv{p) Ji y 



d 



dp p — uo 
(x - yf 



16(p — xY{p — yy 



(3.19) 



(3.20) 



(3.21) 



(3.22) 



(3.23) 
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This enables us to find Wi{p) and we see that it is of the form (^.151) with 



A? = Af^ = 1, (3.24) 

^ 8(i' ^16' ^ ' 

d['^ = i-, Df) = -. (3.25) 

^ 8d' ^16 ^ ' 

Carrying on the iteration process is straightforward. In each step one must calculate 
the right hand side of the loop equation ( p.l2|) . Decomposing the result obtained 



into fractions of the type {p — x) {p — y) " allows one to identify immediately the 



coefficients A'^J'^ and DJ"). 



To calculate Wg{p,p) it is convenient to write the loop insertion operator as 

d ^ dMn d ^ dJj d dx d dy d , , 



dV{p) ^dV{p)dMn ydV{p)dJj dV{p)dx dV{p)dy 

where 

"^^^ ^ ' x)-"-V2 (p _ ^)-3/2 _ + 1/2)<I>("+I)(p) 



dV{p) 




-1)V: 




dV{p) 



(3.27) 



The derivatives dx/dV{p) and dy/dV{p) are given by (|3.21| ) and the function 



was defined in ( |3.11| ). Of course dJn/dV{p) just appears from dMn/ dV{p) by the 
replacements y ^ x and J <-> M. We note that there is no simplification of the 
algorithm in the case of the symmetric potential. We can only put x = —y = \fz 
at the end of the calculation. This complication stems of course from the fact that 
we have to keep the odd coupling constants in the loop insertion operator until all 
differentiations have been performed. Only hereafter they can be put equal to zero. 
The same is not true in the double scaling limit however. We will come back to this 
later. 



By taking a closer look at the loop insertion operator (|3.26| ) and bearing in mind the 
results ( |3.23| ), ( |3.24| ) and ( |3.25| ), it is easy to convince oneself that A^^^ and D^^^ depend 



only on x and y via [x — y] and have the structure shown in equation (|3.1|) . Furthermore 



it appears that Wg{j)\^ . . . ,ps) depends on at most 2 x (3g — 2 + s) parameters. The 
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results for g = 2 obtained with the aid of Mathematica read 



A 



(1) 



(3) 



A 



A 



(5) 



201 



67 Jo 



5 J, 



256 d^Ji^ 
57 



+ 



IIJ2 49 M 2 



49 Jo 



128 rf4 j^-i 32 d3 
^2 



256 Ji 

-^2 



64^5 Ji Ml 
201 



256 Mi^ 
67 M2 



128 d^Ji^Mi 
22 M2 

128 d^JiMi^ 
5M, 



128 c/^Mi^ 
57 

128 Ml 

201 
256^4 Ml 2 " 
67 M2 



+ 



32 d^Mi^ ' 
8 J2 



128dUi^Mi 
3M2 



2 + 



128 # Mr 

'^2 



+ 



49 M 2 



128 ci3 Ji Ml 

5M3 
32rf3Mi3 ' 
5M3 , 67 M2 



256#Mi^ 

J2M2 
64 d^Ji^Mi 



49M22 
' 256 rf2 Mi^ 

J2M2 
128rf2 Ji^Mi 



256 (iMi^ 
201 



+ 



32dMi' 128 d^ Ml' 



256 d^Mi^ 
49 M2 

128 rfM^ 



+ 



15 



5 J2 



128 rf3 Ji Ml 
189 



128 rf2 Ji^Mi 



256 d2 Ml 



2 ' 



105 



256 (i Ml 



2 ' 





= A? 


(M ^ 


J, d — 


-d) , 




= A?) 


(M ^ 


J, d — 






= 


(M ^ 


J, d 


-t^) , 




= 


(M ^ 


^ J, d 




dP 


= 


(M ^ 


^ J, d 





(3.2^ 



The genus two contribution to W{p) is now determined by Eq. ( p.l5| ). 



3.3 The iterative procedure for Fg 

In this section we present an algorithm which allows one to determine Fg, as soon as 
the result for Wg{p) is known. The strategy consists in writing the basis vectors x^^\p) 
and ^("^(p) as derivatives with respect to V{j)). It is easy to verify that the following 
relations hold 

X-<"(rt - (3.29) 
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Combining this with the results ( p.24| ) and ( p.25| ) one immediately finds 

Fi = -— InMi - — In Ji - -Ind (3.33) 
24 24 6 ^ ^ 

which coincides with the expression of Ref. [1^ . 

In the general case things are not quite as simple. The basis vectors can not be 
written as total derivatives. This is of course in accordance with the fact that the 
A and D coefficients now have a more complicated dependence on the potential (cf. 
Section p.2| ). However, a rewriting of the basis vectors allows one to identify relatively 
simply Wg{p) as a total derivative. In the case of x^^Kp) this rewriting reads 

x''"\p) = —' ^ T\n-i-i i ^(i) dy \ ( 1 



Mil 2n-l{r[ ' \ "dV{p)] \x-y^ 

n-l 



2 dM 



2n-l dV{p) 
where "I'l"-' should be written as 



n>2 (3.34) 

k=2 ) 



2n-ltC ^ I ^ 'dV{p)]\x-yl 
, ^ dx 2 dMn-i 

^"•'lyW)- 2n-livU "^2- '^^ 

*i" - (3.36) 

The basis vector x^^Kp) should of course still be written as in ( p. 29 ). The rewriting 
of the \E'(")'s is analogous to that of the x'^")'s. It can be obtained by performing the 
replacements J ^ M and x ^ y in the formulas above. 
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By means of these rewritings we have been able to determine F2. The resuh reads 
119 119 I8IJ2 181 M2 



+ 



7680 J/ rf4 7680 Ml' # 480 Ji' rf^ A80 M^' 
J 2 3M2 



3 Jo 3 Mo llJo' 11 Mo' 



64Ji'Mid3 UJiMi^d? 40 J/ d' 40 M/ci2 
43 M3 43 J3 J2M2 17 



192Mi^d2 ' 192Ji^rf2 ' UJi'Mi'<f 128JiMid4 
21J2^ 29J2J3 35 J4 21 M2^ 
"^160 Ji^d " 128 Ji^ d ^ 384 Ji^ d ~ im Mi^ d 



29 M2 M3 35 M4 
128 Mi^ d ~ 384Mi^c/ 
It is obvious from the formulas above that Fg depends for a non-symmetric potential 
on at most 2 x ["ig — 2) moments. Furthermore, for a symmetric potential, Fg is a sum 
of two identical terms and depends on only at most 3(? — 2 moments. A consequence 
of this doubling for the double scaling limit will be discussed in Section |47^ . 

4 The double scaling limit 

It is easy to determine which terms in the explicit expressions for Fg and Wg{p) de- 
termined in the previous section that contribute in the double scaling limit. However 
it is rather time consuming to determine Fg and Wg{p) away from the double scaling 
limit. In this section we develop an algorithm which gives us directly the result in the 
double scaling limit. Using this algorithm we calculate the correlators and the partition 
function explicitly up to genus four and describe their general structure. 

4.1 Multi-critical points 

Let us consider first the case of the symmetric potential. We hence have x = —y = y/z 
and Jk = {—l)''~^^Mk for all values of k. As is mentioned earlier, the m*^ multi-critical 
point is characterized by the eigenvalue density having (m — 1) extra zeros accumulating 
at both —y/z and y/z, and the condition for being at this point is the vanishing of the 
first (m — 1) moments. For the m*^ multi-critical model the double scaling limit of 
the correlators is obtained by fixing the ratio of any given coupling and, say (72, to its 
critical value and setting 

/ = Zc + an, (4.1) 
z = z,-aA^/™. (4.2) 

The moments then scale as 

J^r^ M^r^ a'^-^ A; e [1, m - 1] . (4.3) 
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Furthermore, it is well known that the genus g {g > 1) contribution to the free 
energy has the following scaling behaviour 

Fg ~ a(2-23)('"+V2) . (4.4) 

Bearing in mind that the structure of Fg is as shown in Eq. ( p.l| ), one finds that the 
following relation must hold 

s I 

J2{m - aj) + J2im - A) - {a + (3)im - 1) > m{2 - 2g) - g - 1 . (4.5) 
j=i i=i 

Since the free energy should look the same for all mult i- critical models, we have 

Nj + NM>2-2g, (4.6) 

^K-l) + E(A-l)<3(7-3. (4.7) 

j=i i=i 



We already know from the analysis in Section that the equality sign must hold 
in ( ^4.6| ). Only terms for which max{s,l) < m and for which the equality sign holds 
also in ([4.7|) will contribute in the double scaling limit. From Eq. (|3.6| ) it follows that 
these terms will have 

7 = ^7-1. (4.8) 

A similar analysis can be carried out for the generating functional. Here it is known 
that the genus g contribution to W{p) has the following scaling behaviour 

Wg{TT, A) ~ a(^-23)(m+l/2)-l (4^g) 

with the exception of Wo{tt,A) which also contains a non-scaling part. From (|3.9|) 
and ( p.lOp it appears that the basis vectors scale as 



Remembering that the structure of the A and D coefficients is as shown in (|3.1|) , one 
finds that the following relation must hold for both A^") and D^"-* 

s I 

E(m - aj) + E(m - pi) - {a + p){m - 1) > m{2 - 2g) + n - g - 1 . (4.11) 

j=i i=i 

again provided max{s,l) < m. Since also the generating functional should look the 
same for all multi-critical models, it follows that A^^^ and D^"^ must satisfy the following 
conditions 

Nj + NM>2~2g, (4.12) 
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EK- - 1) + E(A - 1) < 3(7 - n + 1 . (4.13) 

j=i i=i 

Here we recognize the homogeneity requirement ( pj.lSI ). We know from Section ^.Ij that 
the equahty sign always holds in ( [4.12| ). Only terms for which max{s,l) < m and for 



which the equality sign holds also in ( [4.13 ) will contribute in the double scaling limit. 



We note this means that j = g — 1 for these terms (cf. Eq. ( |3.17| )). 

Let us turn now to the case of the non-symmetric potential, and let us consider the 
^th multi-critical point assuming that the extra zeros accumulate at x. To obtain the 
double scaling limit of the correlators, we first fix as before the ratio between any given 
coupling and the first one to its critical value. Then we scale p and x in the following 
way keeping however y fixed: 

p = Xc + avr, (4.14) 
X = Xc-ah^l'^. (4.15) 

Under these circumstances the J-moments do not scale but behaves as 

Mk ~ aJ^-^ e [1, m - 1] (4.16) 



Using again ([4.4|) we find that the only terms which survive the double scaling limit 



for Fg are those for which 



NM = 2-2g, (4.17) 



s 



^(a,-l) = 3(7-3. (4.18) 
Comparing Eq. ( [4.171 ) with Eq. ( [3.5[) , we see that all the J-dependence disappears 



when the prescription for taking the double scaling limit is as in ( [4. 14) ) and ( [1.15[ ). 
Furthermore, comparing ([4.18[ ) with ( p.6| ) we find that all remaining terms have 7 = 



(7 — 1 as in the symmetric case. The relations ( [4.17| ) and ( [4.18 ) will be exploited in 
Section [^. 



Let us carry out the analysis of Wg{p) also for the non-symmetric case since this 
will give us an additional information about the structure of the D and A coefficients. 
For the basis vectors we have 

^ a— -+1/2 ^ ^(n) _ a-^n . (4.19) 

By exploiting again the known scaling behaviour ( [1.9[ ) of Wg{p), one can just as in the 
symmetric case derive certain homogeneity conditions that the A and D coefficients 
must fulfill. One finds 

s 

: Nm>2- 2(7, E(«i - 1) < 3(7 - n - 1 ; (4.20) 
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4"^: NM>l-2g, ^(a^. _ 1) < 8(7 - 1 . (4.21) 

i=i 

The conditions are different for the A and D coefficients due to the different scahng 
behaviour of the basis vectors. Terms that are important in the double scahng hmit 
are only those for which s < m and for which the equality sign holds in both rela- 
tions in (14.201 ) or both relations in ([4.21|) . Comparing the second relation in (|4.21|) 



with ( [4.13| ) we see that all D-terms vanish in the double scaling limit. Furthermore, 
comparing ( [4.20| ) with (|3.16|) we see that all J dependent A terms disappear. Hence 
in the double scaling limit everything is expressed only in terms of the M's and d, and 
all terms have 7 = (7 — 1. 

As is mentioned earlier, D^"^ can be obtained from A^^"* by performing the re- 
placements X ^ y and J ^ M Hence we also have that the following homogeneity 
requirements should be fulfilled 

Af: Nj>l-2g, ^(/3^. - 1) < 8(7 - 1 ; (4.22) 

Df: Nj>2-2g, j^iP^ - 1) < 3g - n - 1 . (4.23) 

i=i 

These relations could of course be derived by analyzing the scaling behaviour assum- 
ing that the extra zeros accumulate at y. This implies fixing x and replacing (|4.14|) 
and ( ^1.151 ) with 

p = yc-a-K, (4.24) 

y = y, + aA'/^. (4.25) 

In this situation only terms for which I < m and for which the equality sign holds in 
both relations in ( 4.22| ) or both relations in ( [4.23 ) will contribute in the double scaling 



limit. For symmetry reasons, now all A terms plus M-dependent D terms disappear. 
Also the free energy is expressed under these circumstances entirely in terms of the 
J's and d. We note that independently of the details of the prescription for taking the 
double scaling limit, in this limit the multi-loop correlator Wg{pi, . . . ,ps) depends on 
at most 3g — 2 + s moments and the free energy Fg depends on at most 3g — 2 moments. 

4.2 Determining Wg{p) in the d.s.l. 

The homogeneity conditions derived in the previous sections allow us to determine 
whether a given term contributes to the double scaling limit or not. As an example, 
let us consider the case of a non-symmetric potential where the critical behaviour is 
associated with the endpoint x. Using ( [4.20| ) we see that only A'^\ the first term 
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of A^2^ and the first two terms of A^^ of the long hst of comphcated expressions for 



the D and A coefficients of W2{p) given in Section |3.2| survive in the double scaling 
limit. Obviously it would be convenient to have an algorithm which gives as output 
only terms which are relevant for the continuum limit. An algorithm which gives us all 
potentially relevant terms can be obtained by a slight modification of the one described 
in Section By potentially relevant we mean relevant for m multi-critical models 
with m sufficiently large. 

It is easy to convince oneself that terms which are relevant for the double scaling 
limit of Wg{p) can only appear from terms which are relevant for the double scaling 
limit of the right hand side of the loop equation. Let us assume that we know the 
scaling relevant versions of Aq^ and Dq^ for G = 1, . . . ,g and let us assume that we 
want to calculate the double scaling limit of the right hand side of the loop equation 
for Wg^i. First of all we note that all terms in the basis vectors x^^Hp) '^^"'\p) 
show the same scaling behaviour so none of them can be ignored. However we can 
replace all occurrences of {p — y) with (xc — y) = in this limit. It is known from 
earlier analysis fl^ that for a m multi-critical model Wg{j),p) scales as 



Ty^(7r,7r) ~a-2f('"+^/2)-l (4.26) 



Using equation (^4.9| ) it is easy to show that all the products in the sum in (|2.12| ) have 



the same scaling behaviour. Hence, no superfluous terms appear from the sum if we 
start from the double scaled versions of the Wq, G = 1, . . . g. 

However, if we calculate Wg{p,p) by letting the loop insertion operator as written 
in (|3.26|) act on Wg{p), irrelevant terms will appear even if we start from double scaled 



version of Wg{p). By comparing ([4.9|) and ( [4.26|) , we see that only operators in d/dV{p) 



which lower the power of a by m + 3/2 give rise to the relevant terms. Therefore, we 
should discard all other operators. By examining carefully each term in (|3.26|) , one 



flnds that only the following part of the loop insertion operator remains in the double 
scaling limit: 

d _ ^ dMn d dx d . , . , 

'dyW), " ^ dV{p) dMn ^ dy{p}d^ ^ '"""^ ^' ' 

where 



dV{p) 
dx 1 



n + 1/2) {$i"+^)(p) - -^^^^\P)} (d.s.l.) , (4.28) 
^lip) {d.s.l.) (4.29) 



and 



dV{p) Ml 
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We are now in a position to calculate the double scaling limit of the right hand 
side of the loop equation for Wg+i provided we know the double scaled versions of 
Wi{p), . . . ,Wg{p). We see that, since all the y-dependence has disappeared, we do 
not have to perform any decomposition of the result. This is in accordance with the 
outcome of the analysis of the previous section that all D terms disappear in the double 
scaling limit. There it was also found that all J-dependent A terms would vanish. This 
also appears from the formulas above. No J terms ever appear if we do not start out 
with any and we do not. 

The starting point of the iteration procedure is of course as before Wo{p,p) but now 
we should take only the part of it that contributes in the double scaling limit. From 
equation ( |3.23| ) we finJ 

1 1 



wr\p,p) 



16 {p 



x] 



{d.s.l.) 



(4.31) 



There is an interesting feature of the A coefficients. For all genera it holds that A^^-^ = 
A^g^ = 0. This is because the smallest negative power of {p — x) appearing in any 
basis vector is 3/2. By multiplication of two basis vectors or by application of the loop 
insertion operator this power will be lowered by at least 3/2. Hence terms of the type 



x) 



and {p 



X 



-2 



never turn up on the right hand side of the loop equation. 



We have calculated the A coefficients as they look in the double scaling limit for 
g = 2, 3 and 4 following the iterative procedure outlined above. The results for (7 = 2 
and g = 3 read 



A 



(3. 



A 



A^' 



A, 
A 



A 



(3. 



0, 

0, 
49 

256 4 Mi^ " 

-49 M2 
128 4 Ml ^ ' 
105 

256 4 Ml 2 ' 

0, 

0, 

-5355 Ma^ 



5 Ms 



32 dc Ml 



3 ' 



+ 



7995 Ms^ M3 32845 M2 Mg^ + 35588 Ma^ M4 



512rf2Mi9 256^2 Ml® 
21 (680 M3 M4 + 773 M2 M5) 



2048 rf^Mi^ 
1155 Me 



2048 Mi^ 



512 d2 Ml 



5 ' 



^Here and in the following we use the notation that the superscript (5) refers to the case of the 
symmetric potential and the superscript (NS) to the case of the non-symmetric potential (where the 
critical behaviour is associated with the endpoint x). 
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(4) _ 5355 Ms^ 21837 Ma^M: 



^ 2bQdlMi^ 512d2Mi^ 

17545 Mg^ + 37373 Ms M4 10185 M5 



2048rf2Mi'^ 2048rf2Mi^ ' 

7371 Ma^ 69373 M2 M3 17465 M4 



^ ~ 256^2 Ml ^ ^048dp^ ~ 2048dp^ ' 
^(6) _ 64295 Ma^ 30305 M3 

^ ~ 2048dp^ ~ 2048d2M7 ' 
^(7) _ -13013 M2 
^ ~ 512 d2 Mi^ ' 
_ 25025 

^ 2048ci2Mi^ ^ ^ 

We remind the reader that the terms hsted above are only potentially relevant. Whether 
or not they are actually relevant depends on which multi-critical model one wants to 
consider. For a m'th multi-critical model all terms involving M^, k > m vanish in the 
double scaling limit. This is true whether one uses a minimal potential or not. We 
also remind the reader that we assumed that we had a non-symmetric potential and 
that the critical behaviour was associated with the endpoint x. In the case where the 
critical behaviour is associated with the endpoint y all the formulas in this section still 
hold provided dc is replaced with —dc, x with y and M with J. 

Let us turn now to the case of the symmetric potential. In view of the scaling 
relations (^4.1|) and (|4.2| ) it might seem unnatural to work with terms like {p — x) 
and {p — y) and we will see below that there actually exists a way of avoiding this. 
However, for the moment we will analyze the scaling behaviour in the x, y formalism 
by (formally) scaling p to —^/z^ whenever it occurs in a term like [p — y) and to +y/z^ 
whenever it occurs in a term like {p — x). Now let us assume as in the non-symmetric 
case that we know the scaling relevant versions of A^^ and Dq^ for G = 1, . . . ,g and 
that we want to calculate the double scaling limit of the right hand side of the loop 
equation for Wgj^i{p). Of course we still have that all terms in the basis vectors show 
the same scaling behaviour so we still have to keep all of them. However it is easy to 
convince oneself that in the double scaling limit all occurrences of {p — y) in x^"'Hp) 
can be replaced by dc and all occurrences of {p — x) in '^^"'^ can be replaced by —dc- 
Furthermore, it appears that when we calculate the sum of products on the right hand 
side of the loop equation all products that mix x's and \l/'s will be subdominant in 
the double scaling limit. These mixed products should hence be discarded. Finally to 
ensure the survival of only double scaling relevant terms in Wg{p,p), the loop insertion 
operator should be written as 

d d d , , . 

+ TTTTTT (4-33) 



dV{p) dV{p)^ dV{p) 



21 



where d/dV{p)y can be obtained from d/dV{p)x given in equation ( [4.27] ) - ( |4.3CI| ) by 
replacing a; by y, M by J and dc by —dc- 

The starting point for the iterative procedure is of course the double scaled version 
of Wq {p, p) which here reads 

wi'\p,p) = + ^4-34) 

lo [p — xY lb [p — uY 

We see that the loop equation actually decouples into two independent equations. 
Each of these is a double scaled version of the loop equation for the non symmetric 
potential. One corresponds to the case where the critical behaviour is associated with 
the endpoint x, the other to the case where the critical behaviour is associated with the 
endpoint y. The A^"'''s are hence equal to those obtained in the non-symmetric case 
where the critical behaviour is associated with the endpoint x and D^"^ = (— 1)"'A^"^. 
This means that the right hand side of the loop equation of the loop equation can be 
written as 

39-1 1 
the r.h.s. = y d'lA^';^— (4.35) 

which seems to be the natural way of expressing it bearing in mind the scaling rela- 
tions (g3) and (U). 

As a consequence, it becomes also more natural to express the generating functional 

as 

3g-l 

Wf\p)= E4"^X^"nP) (4.36) 

n=l 

where A^") = A^^^d^ and where we now have only one set of basis vectors x^'^\p) which 
fulfill 

{K-2W,{p)}x'-\v) = ^^^^. 
It is easy to show that x^^\p) should be chosen as 



(4.37) 



1 "~ I 

= ^ W:\p) - E X^'\v)Mn-k-,idl-- \ (4.38) 



k=l 



where 



^i"Hp)= (p2_,).^i/2 - (4-39) 

Furthermore, one can easily convince oneself that the right hand side of the loop 
equation can be obtained directly in the form ( [4.35|) if one replaces {p — x)^^ by 
dcip^ — z)^^ in the formulas ( [4.27] ) - ( [4.30| ), i.e. if one carries out the iteration process 
starting from 
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and using the following expression for the loop insertion operator. 



dVsip) V dV{p) dMn dV{p) dz 

where 



-(n + 1/2) {$r+^)(p)rf^ - ^^I'Hp)} 4 (ci.../.) , (4.42) 
"^"^ ^^?kP) {d.s.L). (4.43) 



d\/(p) 



rfy(p) Ml 

This observation allows us t show explicitly that the hermitian and the complex matrix 
model are equivalent in the double scaling limit. 

4.3 Determining Fg in the double scaling limit 

The starting point of the calculation is the double scaled version of Wg{p) obtained 
as described in the previous section — and the strategy is again to rewrite the basis 
vectors in a form which allows one to identify Wg{p) as a total derivative. However this 
time the rewriting in the case of x^'^^ is made with the aid of ( |4.28| ) instead of ( |3.27| ) 
and reads 

X^-Hp) = (:7^^ - j:X^'\p)Mn^k^i] n>2. (4.44) 



Ml [2n-l dV{p) 



k=2 



We do not need any expression for x^^Kp) since the sum in ( [4.44| ) starts at k = 2 
and since we know that in the double scaling limit A^^^ = 0, neither we do need any 
expression for $x(p)- The relevant rewriting of \l/'^"^(p) appears from ([4.44|) when M is 
replaced with J. 

For the non-symmetric potential where the critical behaviour is associated with the 
endpoint x, we immediately find for genus one 

Fi^^^^ = -^InMi {d.s.L). (4.45) 

This expression can of course alternatively be obtained by taking the double scaling 
limit of Eq. (ICTj ). 

Before presenting the explicit results for Fjj^^'^ in genus 2, 3 and 4, let us describe 
the general structure of F^^^^ as it appears from our iterative solution of the loop 
equation. We have 

= E (^1 • • • »s\a, 7);- ^"]^^.;^^"- (d.s.L) g>l (4.46) 
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where the rational numbers denoted by the brackets are identical to those of Eq. ( |3.1| ) 
having the same indices. The a/s obey Eq. ( [4.18[ ), 7 is given by Eq. ( |4.8| ) and 



a = 2g-2 + s (4.47) 

as it follows from Eq. ( |4.17|) . 

Let us consider the g = 2 case as an example of what structures are encoded by 
Eq. ( [4.46| ). FoT g = 2 we have aj G [2, S^f - 2] = [2, 4] and the relation ( |4.18| ) gives 

= 3(7-3 = 3 g = 2, (4.48) 



so that s G [1, 3(7 — 3] = [1,3]. The restriction ([4.48|) admits s = 3 and cti = 02 = as 



2, s = 2 and ai,a2 = 2,3 and s = 1 and ai = 4. The proper powers of Mi in the 
denominator for such terms are given bya = 2(7 — 2 + s = 2 + s. The outcome of the 
iteration process is 

^^NS) ^ -21 29M2M3 _ 35 M4 
^ 160 4 Mi^ 128 4 384 4 ' ^ ' ^ 

We note that all allowed values of s and a^'s actually appear. By taking a closer look 
at the loop insertion operator and the basis vectors it is easy to convince oneself that 
this will be the case for all genera. The expression (|4.49| ) could of course also have 
been obtained by taking the double scaling limit of (p.37|) following the recipe given 
in (|4.17| ) and (|4.18|) . We see that the expression for F2 simplifies drastically in the 



double scahng limit. 

The results for (7 = 3 (where s < 6) and (7 = 4 (where s < 9) read 

^(Tvs) _ 2205 8685 Ma^ Ms 15375 Ma^Ms^ 5565M2^M4 

^ ~ 256d2Mii°~ 256rf2Mi9 ^ bUd^Mi^ ^ 256rf2Mi* 

56O5M2M3M4 72875 Ms^ 3213 Ma^Ms 2515 M3M5 
256rf2Mi^ 21504 rf^Mi^ ~ 256 Mi^ ^ 512rf2Mi^ 
21245 M42 5929 M2 Me 5005 M7 



+ ........ I - o..... '5 (4-50) 



9216 4 Ml'' 1024 4 Mi^ 3072ci2Mi 



and 



(^5) _ 21023793 M2^ 12829887 M2' M3 98342775 M2° M3 



T-^yly'J) ^ 

-^4 — - ,0 , . 1^ + 



^ 10240 rfSMi^^ 1024^3 Ml 4096^3 Mi 

4456305 M2^ M4 16200375 M2^ Ms^ 26413065 M2^ M3 M4 



512rf3Mii^ 1024^3 Ml ^2 1024rf3Mi^2 

12093543 M2^ M5 83895625 M2 Ms^ 68294625 M2^ Ms^ M4 



2048rf3Mi^2 32768 rf3Mi^^ 4096 4 Mi" 
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12367845 M2^M4^ 13024935 Ms'' M3 M5 15411627 Ms^ Mf 



2048c/3Mi" 1024 d3Mi^^ 4096 rf^ Mi" 

32418925 Ms^ M4 17562825 M2 M3 578655 MaMg^M 



+ ^^^r^TTTln + 



24576 dl Ml ^° 4096 dl Mi ^° 1 28 d^ Mi ^° 

10050831 M2^ M4 M5 5472621 Ma^ Ms Me 44207163 Ms^ M7 



2048 dl Ml 1024 df Mi 20480 Mi ' 

5 



1511055 M2M52 7503125 M4^ 2642325 M3M4M 



2048rf3Mi9 36864 rf^Mi^ 2048 d^Mi^ 

11532675 Mg^ Me 6242775 M2 M4 Me 6968247 M2 M3 M7 



16384 d^Mi^ 4096 d^Mi^ 4096c/3Mr 

4297293 Ma^ Mg 12677665 M2 Mg 8437275 M5 Mg 



4096 d3Mi^ 32768 d3Mi« 32768 rf^Mi^ 

8913905 M4M7 10156575 M3 Ms 8083075 Mio 
^ 32768 rf^Mi^ ^ 32768^3 Ml ^ 98304 rf^ Mi ^ ' ^' ' 

Needless to say that the results for the case where the critical behaviour is associated 
with the endpoint y can be obtained from these by performing the replacements M — 
J 5 dc ^ dc. 

It is obvious from the discussion in Section ^?1] and formula ( [1.44| ) that indepen- 



dently of the details of the prescription for taking the double scaling limit, Fg depends 
in this limit on at most 3g — 2 moments. Again we stress that the terms listed above 
are only potentially relevant. For a m'th multi-critical model all terms involving M^, 
k > m vanish too. This is true whether one uses a minimal potential or not. It is in- 
teresting to note that for the symmetric potential we get a sum of two identical terms 
and hence 

Ff ) = 2F^(^^) {d.s.L) (4.52) 
which is a well known property of the double scaling limit for the hermitian matrix 



model [0. It can be traced back to the fact that the loop equation for a symmetric 



potential decouples into two independent equations. Equation (|4.52|) is the reason why 



the partition function of 2D quantum gravity is defined as the square root of the double 
scaling limit of the partition function ( p.l|) in the case of a symmetric potential. There 
is no such square root in the non-symmetric case and the proper continuum partition 
function can be obtained just as the double scaling limit of (|2.1|) . 

Of course we have determined up till now only the coefficients, Fg, of the genus 
expansion of the free energy (Cf. Eq ( |2.11D ). For a m'th multi-critical model the 
relevant expansion parameter in the double scaling limit is iVa™^^/^. If we define 
moments fik by (Cf. Eq ([4.3|) ) 

Mfe = a'"-Vfc, ke[l,m] (4.53) 
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we get by replacing Mk with /i^ for A; G [l,Tn] and setting Mk equal to zero for k > m 
in the formulas above exactly the coefficient of the expansion in the double scaling 
parameter. A similar statement holds for the results of the previous section. 



4.4 A remark on the complex matrix model 

The complex matrix model is defined by the partition function 

ZC^^N^pc^f d(j)U(j)exp{-NVc{(j)^(f))) (4.54) 

JNxN 

where the integration is over complex N x N matrices and 



oo 



Its generating functional is defined by 

1 

N 



^c(0V) = E-Tr (0W- (4.55) 

j=i J 



1 oo 

W'^ip) = ^E(Tr (</>V)')y'+' (4.56) 

k=0 



and the s-loop correlator {s > 2) by 



iy^(pi,...,p,) = iV-2 E (Tr(0V)^^..Tr (0^)^=0 Un./p^+^..p; 

fci,...,fcs = l 



(4.57) 



PF^(pi, ...,ps)= N^-' ( Tr ^ . . . Tr ) . (4.5^ 



so that 

\ Pi -</)'(/> Pt 
In Ref. [|l^ an iterative procedure which enables one to calculate W^{pi, . . . ,Ps) 
for any genus g and any s starting from W^{p) was developed. This procedure is much 
the same as for the hermitian matrix model. It is based on the loop equation, the 
loop insertion operator and a suitable change of variables. In Ref. [|19]| the studies were 
carried out only away from the double scaling limit. In the Appendix we show how 
the iterative scheme for the complex matrix model can be modified to give directly 
the result in the double scaling limit. The modified scheme provides us with an exact 
proof that the hermitian and the complex matrix model are equivalent in the double 
scaling limit!. We find that 

W^ipu...,Ps) = -^,W^'\pi,...,Ps) (d.s.L), (4.59) 

= ^Ff) id.s.l) (4.60) 

^The fact that the hermitian and complex matrix models belong to the same universality class and 
that the correlators therefore coincide in the double scaling limit docs not contradict Ref. ||2^ where it 
was shown that the discrete Virasoro operators for the complex matrix model do not have continuum 
limit. 
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independently of the type of critical behaviour considered. In particular we see that 
the partition function of 2D quantum gravity can also be obtained as the square root of 
the double scaling limit of the partition function ( [4.54 ) provided the integration is over 



N/2 X N/2 matrices in stead of over N x N matrices. This is of course in accordance 
with the fact that the number of independent components of a N/2 x N/2 complex 
matrix is the same as that of a x hermitian matrix. 

5 Relation to moduli space 

The equivalence between the double scaling limit of the hermitian matrix model and 



the Kontsevich model has been proven in various ways in Refs. [^6|, We present 
in this section another proof directly in terms of the moments introduced above. For 
this purpose we describe a limiting procedure which allows us to obtain the Kontse- 
vich model directly from the hermitian matrix model. We discuss also the connection 
between the hermitian matrix model and a certain discretization of moduli space. 

5.1 The limiting procedure 



As we have seen in Section ^.3| , the double scaling limit of the free energy is given by the 
expansion ( [4.46|) which encodes in a remarkable way all multi-critical behaviour of the 



hermitian matrix model. Since most of the terms originally present in the expression 
for Fg vanish in the double scaling limit the question arises whether it is possible to 
formulate a matrix model which gives solely the coefficients which are relevant for the 
double scahng limit The answer to this question is — yes — and the proper model is 
known as the Kontsevich model. 

The Kontsevich model is defined by the integral ^ 

/„,„<iA'exp{ntr {-'^MX')} 

there M=diag(mi, . . . ,m„) (at first) is a positive definite matrix and the integration 
is over hermitian n x n matrices. The genus expansion of the free energy is given by 



CO 



In Zk[M, n] = J2 n^-^^F^°''' . (5.2) 

3=0 

The goal of this section is to show that the Kontsevich model can be obtained from 
the hermitian matrix model by a certain limiting procedure which looks conceptually 
different from the double scaling limit. We shall see however in Section |5]^ that this 
limiting procedure leads to a result which is identical to the one obtained in the double 
scaling limit. 
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In order to formulate the limiting procedure let us recall some identities. Consider 
the hermitian matrix model as defined in Eq. ( |2.1| ) with N = —^riE We have the 



following identity |18 



Z[g,N{0] = e-'^'^^'/'ZplA, 



n\ 



iV(0 = -^n 



(5.3) 



where the partition function Zp[A,n] is associated with the Penner model |]29|, in the 
external field A =diag(Ai, . . . , A„): 



Zp[A,n] = f dXexp [ntr (aX - - ^logx) 

Jnxn L V 2 / 



(5.4) 



The identification is valid provided the set of coupling constants {g} = {go, gi, . . .} and 
the matrix A are related by the Miwa transformation 



gk 



-tr A"'^ 
n 



k2 



k> 1, 



9o 



-tr log A ^ 

n 



(5.5) 



In order to show how the partition function ( p.l| ) can be obtained from ( p.4| ), we 
note that Zp[A, n] can be reformulated as the so-called Kontsevich-Penner model 



ZKp[i],n] 



/nxnC?^exp{-^n tr 


'^r]Xr]X + {\og{l-X)+X) 


} 


InxndXexp[-^ntT 


'^r]Xr]X - |X2 


} 



provided 



(5.6) 



(5.7) 



We note in addition that the model (|5.6|) possesses a remarkable connection to a dis- 
cretization of moduli space which will be discussed in Section |5.4| . Here we will be 
interested in the connection to the usual continuum moduli space, i.e. to the Kont- 
sevich model which precisely has the interpretation as the generating functional for 
intersection indices on moduli space (Cf. to Section 5^). In order to highlight this 
connection we choose a specific parametrization of rj and ^: 



V 



^ 2e3 



(5. 



where m is to be identified with the matrix, M which appears in the definition (|5.1|) 
of the Kontsevich model, and where e is presently just an expansion parameter which 



^We should point out a subtlety related to the representations (5.4) and ( ^.6[ ). We have used a 
constant ^ in these equations which has the opposite sign of that conventionally used. The "analytic" 
continuation of ^ is needed in order to make contact with the Kontsevich model. We will not to try 
to "justify" it, since the idea is to provide some heuristic explanation of the observed identity of the 
expansion coefficients of the hermitian matrix model and the Kontsevich model. 

®See also Ref. 28 where the proof was extended to finite n. 
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we are going to take to zero, but which in Section will be given the interpretation 
of the step of discretization of moduli space. 

It is now trivial to check that we in the limit e reproduce (SJ.) from (|5.6|) 
after rescaling X eX. In this way we have managed to move continuously from 
Z[g,N{C,)] to ZK[iTL,n]. Note that there is no need to take the limit n — ^ oo during 
these steps, but the size N{C,) = —C,n of the original hermitian matrix goes to infinity 
as e — s> 0. 

It is instructive to illustrate the proof just given by taking the e —>■ limit of the 
explicit formulas governing the hermitian matrix model. The basic equations (|2.14|) 
which in the original matrix model determined the endpoints x and y in terms of the 
coupling constants Qi have an equivalent formulation in terms of the eigenvalues Xj of 
A when we rewrite the matrix model in terms of the Penner model as in (|5.3| ). They 
are obtained by inserting in ( p^.l4 ) the following relation between the eigenvalues Xj 
and the matrix model potential, V: 

v-'M^ig^— (5.9) 

This relation follows directly from (|5.5|) . By the use of ( |5.9|) the equation ( p.l4| ) reads 

1 " 1 x + y , , 

-Y , — = 0, 5.10 



^J{Xi-x){Xi-y) 



«S^(Ai-i)(A,-j,) 8 

Notice that the normalization factor "2" in ( |2.14| ) has been changed to —2^, since the 
size of the matrix is = — ^n. 

Let us now apply the limiting procedure dictated by (|5.8| ) to Eqs. ( ^.lOj ) and ( |5.11 

V2 



A 



y^(e - + e--) = ^ + y|m2 + 0{e'/'). (5.12) 
In order to solve equations ( |5.1CI| ), ( |5.11| ) in terms of x and y we choose 

l/=-||, x = ^^+V2e'/\o + .-. . (5.13) 
From Eq. ( |5.10| ) we have the equation 

It : ' =Uo (5.14) 

which appears multiplied by y^. After the substitution into Eq. ( |5.11| ), we get that the 
leading term of order coming from 2^ and (x — ?/)^/8 cancels. It is crucial for this 
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cancellation that ,^ > 0. The next order terms which are proportional to 1/e combine 
again into Eq. ( |5.14| ), which we recognize as nothing but the stationary condition 
appearing in the Kontsevich model. 

Let us now demonstrate how the 1-loop average of the hermitian matrix model is 
related to the one of the Kontsevich model. We start from genus zero inserting Eq. ( |5.9| ) 
into (|2.13| ) which gives |jl8 



Woip) 



ip-x){p-y)-p+-J2 



{p-x){p-y) 



- 1 



(5.15) 



n , p-\ \^^X^-x){\,-y) 

Eq. ( ^.10|) ensures the vanishing of the 0{1) term of the 1/p expansion of the right 
hand side, while Eq. ( p. Ill ) implies that the 0{p~^) term equals —^/p- Contrary the 
representation ( |2.21|) can be obtained from (|5.15|) by expanding in 1/Aj and using 
Eqs. (|3D and (p|). 

To obtain an analogue of Eq. ( p.l2| ) for p, let us point out that the relation ( |5.5| ) 
implies 



d_ 



1 d 



E 



d 



n 



oo 



(5.16) 



dgo t^i +^ dgk 

which requires n — oo in order for (yf^'s to be independent variables. Comparing with 
Eq. (p^) one gets 

n OO . (5-17) 



dV{p) 



_d_ 



Therefore we substitute 



P 



72 

^3/2 



+ 



/f- + 0(.^/^) 



(5.1^ 



where tt is to be understood as the momentum variable for the Kontsevich model. 
Substituting (|53^) , ([5J3| ) and (|]T|) into Eq. ( ^J5D , one gets 



w,{p) + 1 



2n Y P - \ 



2n Y 



ml 



(5.19) 



where 



konti 



d 



1 
2 



2 lnZo^°"* 

m-'=-K 



[M, n] 



Vvr - 2mo - + - XI 



n J IT — TTlj 



2Un 



rrin 



(5.20) 



We see that the subtraction needed to make contact with the Kontsevich model is 
exactly the same as the usual double scaling subtraction. For the two-loop correlator 
we have 



Wo{p,q) + - 



1 1 



2 (p - g)2 



2(vri - 7r2)2^7ri7r2 



(5.21) 
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We recognize the term l/2(p — qf as being the non universal part of the two-loop 
correlator. Hence the renormalization needed in this case is also the usual one of the 
double scaling limit. For p = q we find 



'2 

Wo{p,p) 



^"'"^^16^J '^^•^^^ 

The appearance of the term I/IGtt^ when one wants to make contact with continuum 
physics is well known from the study of the double scaling limit of the Virasoro con- 
straints for the hermitian matrix model [^. It is evident from equation (|5.21| ) that 
for the higher loop correlators it holds that 

Woipi, ...,Ps)^ {^^y Wt'^'i^u • • • , vr.), s > 2 (5.23) 

This multiplicative renormalization is also in compliance with the renormalization of 
the double scaling limit. The origin of the factors (2/e)^''^ is easy to understand bearing 
in mind the relations (|5.12| ) and (|5.17|) . In the next section we show that the limiting 
procedure described above actually leaves us with the same terms as the double scahng 
limit does. It is worth noting that this means that the Kontsevich correlators encode 
information about continuum correlators at all multi-critical points. This information is 
even easily accessible. To obtain the continuum correlators one rewrites first Eq. ( |5.14| ) 
as the genus zero string equation 

E ^^^^^^^0 = (5.24) 



where 



tk = -tT M-^'^-i _ Ski k>0 (5.25) 

n 



are continuum coupling constants which all vanish at the m^^ multi-critical point ex- 
cept for tjn and to with the latter playing the role of the cosmological constant. The 
continuum correlators can hence be found by expanding the Kontsevich correlators in 
powers of l/frij. Equation (|5.14| ) ensures the vanishing of the l/-\/vr term in the the 
expansion of the right hand side of ( ^.20| ). 

We are now ready to present one more proof of the fact that by carrying our limiting 
procedure for the hermitian matrix model we recover the Kontsevich model. For this 
aim let us consider the Schwinger-Dyson equation of the model ( |5.4D which reads 

1 dW{\, 



11? dXj 



where 



W(A.)^,y(A.) + | + i^i:3^ (5.27) 



31 



Eq. ( |5.26| ) can be obtained by inserting Eq. ( |0| ) into the loop equation ( |2.12| ). 

Substituting Eqs. (|5.12| ), ( [5.20| ) and (|5.22| ) and assuming multiplicative renormal- 
izability of the higher genera contributions to the correlators (as in the double scaling 
limit) one gets from equation ( 5.2(j| ) the following equation 

9 1 y{/kont( 



1 aiy'=°"*('m?l 



where 



dm'j 



n '-r'. m?, — mf 



4 



— H V — 

2 2n^.mi 




This is exactly the Schwinger- Dyson equation of the Kontsevich model. 



(5.2^ 



(5.29) 



5.2 The equivalence of the d.s.l. with Kontsevich model 

The equivalence of the Kontsevich model and the double scaling limit of the hermitian 
matrix model in genus zero was demonstrated in Ref. |^] by comparing explicit so- 
lutions. In order to compare next orders of the genus expansions, let us consider the 
following expression for conjectured in Ref. 0] for the Kontsevich model: 

mt ^ai ' ' ' In 

'a 



-iKont 



^ (ai . . . as\a) 



9>l 



where the moments /i,'s are defined b 



(5.30) 



4(M) 



n 



E 



{ [m] - 2uo)^+i/2 



A; > 



(5.31) 



and uo{M) is determined by the equation uq = Io{uq, M), i.e. by Eq. ( p.l4| ). When we 



compare the coefficients calculated in Ref. JT^ with the i^.^"-'^^™- calculated in the 



, herm 



double scaling limit in Section |4^, we see that they are identical. In this section we 
show that this is true to all orders in the genus. 

For any value of e in (|5.8| ) we have the expansion ( p.l|) of Fg of the hermitian 
matrix model in terms of M^'s and J^'s. We are now going to show that the terms 
which survive in the limit e ^ are the same as those that survive in the double 
scaling limit ( [4.46| ), and that the M^'s are related directly to the J^'s of the Kontsevich 
model defined by ( 5.31 ) and ( ^.14 ). 

Let us now discuss the scaling behaviour of the moments Mk and Jk- By substituting 
]9D into Eqs. (|2.171) and (|2.18| ), we obtain for the moments M^, J^: 

1 



n 



(A, -x)^'+i/2(A^. -y)l/2 



Jkl 



k> 1 



(5.32) 



^This definition differs by a factor —{2k — 1)11 from the Ik's used in Ref. 0. 
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From this representation and Eqs. ( p.l2| ), ( p.l3| ) it is easy to infer for the moments 
and Jfc and for the parameter d the following scaling rules: 

Ju ^ _2-(3fc/2+l)^(3/c+l)/2j^ ^ ^^^^ (5^34) 

M, ^-2(^-i)/V(^-^)/2(4-5,0, (5.35) 
d ^ 23/2^-3/2. (5.36) 



where /^'s are the the Kontsevich moments defined by (|5.31|) , (|5.14|) . 

Let us turn to the general expression ( p.l| ) for Fg. Then the power in e for a specific 
term 

is 



Clearly, due to the formulas (p^ ) and (|3.8|) , we find that [e] > 0. [e] = is possible if 
and only if / = 0, 7 = (7 — 1 and 

^(«,— l) = 3(7-3. (5.39) 

This means that we have precisely the same terms surviving as in the generic double 
scaling limit ( [4.46| ). By taking a closer look at the loop insertion operator ( |3.26| ) one 



finds that such terms always have /5 = and hence Nm = 2 — 2(yf. We therefore 
conclude that starting from the generic hermitian matrix model and letting e — we 
will reproduce ( |5.30D , ( p.31| ) and ( p. 141 ) the coefficients in front of the Ik terms being 
those inherited from the hermitian matrix model. Hence we have = (■)^°"* since 

we saw explicitly that the partition function Z[g,N{S,)] is identical to the partition 
function of the Kontsevich model in the limit e — 0, ^ = l/2e^. 

Thus the coefficients which survive in the double scaling limit have an in- 

terpretation as intersection indices on moduli space, since this is the interpretation of 
the coefficients of the Kontsevich model. In Section |5.4| we are going to show that the 



additional coefficients which appear away from the double scaling limit can also be 
given a geometric interpretation on moduli space. 

5.3 Intersection indices on moduli space 

It has been known for some time that matrix models might be very useful in providing 
explicit realizations of the moduli spaces 7Wp,s of Riemann surfaces of genus g and s 



33 



marked points Xi, . . . ,Xs- The first example is the Penner model [EUl, which allowed 



a relatively simple calculation of the virtual Euler characteristics of these spaces. The 
next example is the Kontsevich model, which allows a calculation of intersection indices 

on Mg^s- 

It is not our intention here to define or discuss these concepts in any detail. Let 
us only note, to fix the notation, that there exist s natural line bundles Ci, i = 1, s 
on Mg^s- The fiber of £j at a point S G Mg^s is the cotangent space to the point x,, 
viewed as a point on the surface S. The line bundles have first Chern classes ci(£j), 
which can be represented by the curvature 2-form of an arbitrary U{1) connection on 
Ci- If we have non-negative integers such that 

« 1 

^ ctj = -real dim A4g,s = 3g — 3 + s (5.40) 

1=1 

it is possible to form the integral [0] 

(ra,---r«.)= /_ ci{C,r---c^{Csr (5-41) 

where products and powers are exterior products and powers and ^Ag^s refers to a 
suitable compactification (the so-called Deligne-Mumford compactification) of M.g,s- 

The numbers (tq,^ ■ ■ ■ Tq^, ) are topological invariants and one of the achievements of 
Kontsevich was that he found a more manageable representation by using a combi- 
natorial decomposition of M.g,s inherited from the physicists "fat-graph" expansion of 
the hermitian matrix integrals. If we have a double-line graph with s faces (dual to the 
s punctures), we can assign a perimeter pf = J2eh to each face, /g being the length of 
the edge e and the sum being over edges which constitute the face /. The set of such 
double-line graphs with assigned Vs is a combinatorial model of moduli space. 

Cells G have dimensions over the real numbers obtained by counting the number of 
independent l^s for fixed p/s. The result is that 

real dim G < 2{3g - 3 + s) (5.42) 

where the equality sign is valid if and only if all vertices of the graph G are trivalent. 

For each face ft one can introduce the 2-form Ui, which can be considered as the 
combinatorial version of ci(£i) related to the puncture xf. 

co^ = T.d(—]^4—] (5-43) 

a<b \Ph) \PfJ 

where la are the lengths of the edges of fi, which are assumed to be cyclic ordered. For 
fixed p/s we can now write the formula for intersection indices as 

{rar--raj = (5.44) 
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provided the integral is over cells G of maximal dimension, i.e. provided all vertices 
are trivalent. Let us define denote the complex dimension of Aig^s by dg^s, i-G. 



dg,s = 3(7 — 3 + s , (5.45) 

and introduce the 2-form 

n = j2pluj,. (5.46) 

1=1 

It acts as a generating function for intersection indices: 

„ s „2«i 

y^'^'^= E {^^.■■■rO..)gdg,\\{'-^ (5.47) 

ai + ...+ois=dg^s i=l 

By a Laplace transformation we can trade the pj's for mj's by 

dp, e-P^™^ p'"' = (2a,)! m"'"'-^ (5.48) 
and one of the main results of Kontsevich is that 

where the summation on the right hand side is over all double-line, trivalent, connected 
graphs G of genus g and with s faces. The quantity denoted by V is the number of 
vertices in G, e denote an edge in G and / and /' the faces sharing a given double-line 
e. Finally Aut{G) is the order of symmetry group of G. For future use we note that 
that the last factor can be written as 

n ---- = / n die e-'^(-/+-/') (5.50) 
^^^mj + mf, Jo 

When we compare the right hand side of Eq. ( p. 49 ) with the definition of the 
Kontsevich integral (|5.1| ) we see immediately that ( ^.49[) represents the sum of the 
connected diagrams of genus g and s faces which would be generated by a perturbative 
expansion with respect to the gaussian part of the action. In the next section we shall 
extend this result to a more general interaction. 



5.4 The discretized moduli space 

As we shall see the results obtained in this article away from the critical regime is 
related to a representation of moduli space which incorporates explicitly the boundary 
components of the Deligne-Mumford compactification (or reduction) procedure. The 
basic ingredient in the new construction is a discretization of A4g°Y^^. The length Zg 
of an edge is assumed the belong to eZ^, where e is an expansion parameter, which 
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eventually is going to be associated with the one in the Section f).'2\ . This means that 
the Laplace transform in pi performed in ( |5.48| ) is going to be replaced by a sum. The 



effect of this summation is heuristically described by referring to the representation 



pOO I 

/ n die e-''^^™^ ^ e^d^.. -Q ^ ^-en^mf+m'^) ^ ^3d,,. (5.51) 

eeG eeGne=l fo) ™J - 1 

where the power of e is valid for connected trivalent graphs. In the last product we 

recognize the propagator for the Kontsevich-Penner model ( p. 6] ) which involves the 
factor 

^ (5.52) 



ViVj - 1 

provided we make the assignment rji = exp(emj). This assignment is nothing but the 
one already made in ( |5.8|) , and the e assignment in (|5.8|) is also consistent with the 
e factor in ( p.51[ ). In order to fit our notations to the ones by Kontsevich it will be 



sometimes conveninent to use = 2^ = e ^. In fact the perturbative expansion of the 
Kontsevich-Penner model (15.61) can be written as 



oo oo 9_9 _ o— V 9 

1. Z,,lri, „] = E E n'--^' (?) ' ' E n (5.53) 

g=Qs=Q G ^^^V-^) {ij)VtV] J- 

where the notation is as in ( ^.49| ), except that Y1,g denotes the summation over all 
connected double-line graphs of genus g and s faces, not only the trivalent graphs as 
in (|5.49|) . In this way we see that the Kontsevich-Penner model for finite e provides 



a generalization of the Kontsevich model in the sense that it in the limit e ^ just 
reduces to the Kontsevich model, while it for e > describes a discretized version of 
moduli space, which however allows us to access the boundary in more detail. 

By use of this formalism it can be shown that the Kontsevich-Penner model allows 
an expansion, related to the Deligne-Mumford reduction in moduli space, which can be 
viewed as a generalization of the corresponding expansion of the Kontsevich model ||2^ . 
For the coefficient of n^^^d^'^-'ig-s -^^ expansion ( ^.53|) of In Z^p one has 



'6+3s 



E 



reductions 



d, 



E 

Y^aa=dj 



Tai ■ ■ ■ T^af.TQ ■ ■ ■ Tq 



(/.)! 

9, {f, + k,)\ 



n 7 



t}^\edmk) nli(e^"^^-l) 



(5.54) 



symmetrized 
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As opposed to what was the case for the Kontsevich model one should now take into 
account all reductions of the original Riemann surface. The notation is as follows. 
For a given reduction, p denotes the number of connected components of the reduced 
Riemann surface. These components are labeled by the index j. The j'th component 
has genus gj, fj original punctures and kj punctures resulting from the reduction 
procedure. Note that there is no variables corresponding to these new insertions, 
that leads to appearing of additional symmetrical factors + \ in ( ^.54|) . The 

corresponding moduli space has complex dimension dj = Sgj — S + fj + kj. The quantity 
Vp is the power of reduction, defined by J2^=i dj = d ~ Vp. The symmetrization ensures 
that Fg is invariant under the replacements — > (— 1)'^+^ as it must be. 

The expression ( [5.54| ) resembles in many details the original answer by Kontsevich, 
the only two differences being that firstly the variables of expansion, which in the 
Kontsevich case appear in the combination — (2aj — 1)!! tr m~^"*~^, in ( p.54| ) appear in 
a more involved expression. Secondly the complete answer ( p.54| ) contains information 
about the reductions of the Riemann surface. 



Since, as shown explicitly in Section the Kontsevich-Penner model can be 
mapped onto the general hermitian matrix model the information about this reduction 
is encoded in our expansion of the hermitian matrix model in terms of the moments 
Mfc, Jk, and by exploring it we may hope to gain some information about the boundary 
of the moduli space. In order to do this, we need to expand moments M^, Jk and d in 
terms of the somewhat unusual expansion parameters standing on the right hand side of 
Eq. ( |5.54| ). As the first step we need to resolve the constraint equations ( |2.14| ) in terms 
of these parameters. As the first approximation we have (in terms of eigenvalues): 

~-;^ng(l-e-0^^+ni+e^-) ^^-^^^ 

^ n g£mi{fc+l) 

^~Vf^^S(l-e--0(l+e^-0''+' ^^'^^^ 
d ^w^/4--iV- A -\. (5.57) 

Let us now check the validity of the expansion ( p.54| ), for the simplest non-trivial 
example, the moduli space Fi^i (a torus with one marked point). Since s = 1, we should 
extract only single trace terms from the formal expansion of the answer in product of 
traces. Using the explicit answer for Fi, we find that 



■— IndViMi 
24 



tr 



24^ [ (e*^"^ + l)(e^'" - 1) 

.(5.5^ 



■gem ^ X)3(l - e^™) (e^'" + l)(l-e 



£m\'i 
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Doing this sum we reconstruct the answer (which can be also obtained independently 
using diagram technique): 



1 3 e 2^"^ - 1 
^1,1 = -T7T7tr 



tr 



12^ (e2^™-l)3 
92 1 



24^ £2 e ^"^ - 1 e ^"^ - 1 



symmetrized 



(5.59) 



where we again substitute ^ = 2^ = e ^ in order to compare with the Kontsevich model 
answers. That is, the form ( 5.54 ) is explicitly reproduced. The first coefficient, ^, is 



just the Kontsevich index, and the second term corresponds to the unique reduction 
of the torus, i.e. the sphere with an additional factor 1/3! due to reduction. 

This rather simple, but instructive example shows that the hermitian matrix model 
might be able to offer some new insight in the fine structure of moduli space. 

6 Discussion 

We would like to mention at the end one more approach to the hermitian matrix model 
which is based on its interpretation as an integrable system. The appropriate hierarchy 



for the hermitian matrix model is the Toda-chain one [^, ^] and its partition function 



can be represented as the corresponding r-function. This r-functin obeys |Q a set of 



discrete Virasoro constraints quite similarly to the continuum partition function |34| 



The relation between these two sets of the Virasoro constraints was studied in Ref. 25 



As is discussed in Ref. [^, these Virasoro constraints for the hermitian matrix 



model are nothing but the loop equation ( p.l2|) . Therefore, the explicit solutions found 



in the present paper are simultaneously solutions of the Virasoro constraints. We 



do not refer, however, to the fact that the partition function (2J.) is the r-function. 
In the double scaling limit this relation gives a powerful computational method to 
obtain higher orders of the genus expansion based on the Gelfand-Dikii technique. Its 
analogue for the hermitian matrix model away from the double scaling limit is yet 
missing. We think that our explicit calculations prompt that such a technique should 
exist. 
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Appendix A The d.s.l. of complex matrix model 



The partition function of the complex matrix model and the definition of its correlators 



were given in Section O. In the same way as for the hermitian one matrix model, 
one can obtain the correlators from the free energy by application of the loop insertion 
operator, 

l^^(pi,...,p,) = / — ^ (A.l) 

dVciPs) dVc{Ps-i) dVcipi) 

where 

d °° i d 

In Ref. |]T9[ an iterative scheme for calculating W^{p) for any genus g starting from 
Wq{p) was developed. The idea was very much the same as described in Section |3^ . 
The basic ingredient was the loop equation for the complex matrix model: 



which after insertion of the genus expansion of the correlators reads 

d 

dVcip) 



{kc - 2<(p)} W^\p) = E W^ip) W^_^,{p) + ;^^W^-iip) ■ (A.4) 



The operator Kc is given by 

KcHp) ^ I ^ (A.5) 

Jp Am p^ — uj^ 

where Vciyj) = J2j9j^'^''/j and P is a path which encloses the singularities of Wc{p)- 
It was assumed that these consisted of only one cut [— V^? on the real axis which 
in the eigenvalue picture corresponds to assuming that the eigenvalue density for 
has support only on the interval [0, z]. 

To characterize the matrix model potential, instead of the couplings gj, the following 
moments were introduced 

4 - k>0. (A.7) 



P 



4m (w;2 + c)'=+V2 



For the complex matrix model the m^^ multi-critical point is characterized by the 
eigenvalue density for 0^ having a zero of order (m— 1) at the endpoint z of its support. 
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It is easily shown that the condition for being at this point is Ji = /2 = . . . = Im-i 

9 



and Im ^ 0. In Ref. it was proven by induction that W^{p) can be written as 



39— 1 g 

n=l m=l 

The basis vectors X(")(p) and Y^'^\p) fulfill 

{kc-2W^{p)}x^-\p) = (A.9) 

{i^c-2<(p)}r(")(p) = ^ (A.IO) 

and are given by 

= , (A.ii) 



Y^-\p) = -^\n^-\p)-j:Y(''\p)N^.,\ (A.12) 



fe=i 



where 



The iteration process for determining W^{p) is quite similar to the hermitian matrix 
model. To start one needs an expression for W^{p,p). Using a trick similar to ( |3.22|) , 
one finds 



In the general step one makes use of the following rewriting of the loop insertion 
operator 

d ^ dlj d ^ dNj d dz d 



where 



dVcip) ^dVc{p)dh ydVc{p)dN, dVc{p) dz 
dh 



dVcip) 
dNj _ 1 



_ (, + 1/2) - :^$W(j9)| 2 , (A.17) 



{2j + l)n^^^'\p) 



dVc{p) 2 \ ' ' '"^^ p^^{p^ + c)^/'^j 

did) - i-r^^^- ^^-''^ 
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Calculating the right hand side of the loop equation one gets an expression involving 
fractions of the type — z)~"^. Decomposing these fractions into fractions of the 

type (p^ — z)~^ allows one to identify immediately the coefficients A^J^^ and T)'^"'\ 
Let us now modify the iterative procedure so that it gives us directly the result 
in the double scaling limit. The double scaling limit of the correlators for the m^^ 
multi-critical model is obtained by fixing the ratio of any given coupling and the first 
one to its critical value and scaling p"^ and z as in ([4.1|) and ( [4 .21) . It is easy to convince 
oneself that, just as it was the case for the hermitian one matrix model, adjusting the 
iterative procedure to give only the double scaling relevant results amounts to starting 
from the double scaled version of Wq{p,p) and to discarding all operators in d/dVc{p) 
which do not lower the power of a by m + 3/2. One finds 

<(P^P)= 16(/_,)2 ^d.s.l.) (A.20) 

and the loop insertion operator reduces to 

d dli d dz d / , , X / A N 

T.:JT7:7-,l^ + :JT7:7-,ir id.s.l.) (A.21) 



dVcip) Y dVcip) dh dVcip) dz 



where 



dli ,. 



1/2) - (d.s.l) , (A.22) 



dVcip) 

= Y<I>''\P) id.s.l.) . (A.23) 
dVcip) h 



Let us compare the formulas ( |A.2C ) - ( A.23 ) with the corresponding formulas for 



the double scaling limit of the symmetric hermitian one matrix model ( |4.4CI| ) - ( 4.43 ). 
Noticing that d^ = Azc and that in the double scaling limit 



Mfc = <-^4, (A.24) 

we see that in this limit W^{p,p) = 1/4 ■ wi^\p,p), d/dVc{p) = 1/4 ■ d/dVs{p) 
and X^"'\p) = x^'^\p). Since the loop equation looks the same for the complex and 
hermitian matrix models (cf. Eqs. ( |2.15|) and (|A.4| )), we thus have 



W^f (Pi, ■■■,Ps) = ^^,W^'\pu ■■■,Ps) (d.s.l.) . (A.25) 

Furthermore, since the generating functional is related to the free energy in the same 
way for the two models, we also have 

= ^Ff ) (d.s.l.) . (A.26) 
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Hence, the hermitian and complex one matrix models are equivalent in the double 
scaling limit. However, away from the double scaling, we clearly have to distinct 
models. We draw the attention of the reader to the fact that no assumptions about 
the type of critical behaviour entered the arguments of this appendix. 
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